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Introduction 

In  March  1986,  NAVOCEANSYSCEN  tasked  ISA  to  prepare  a  comprehensive  test 
suite  to  systematically  test  and  evaluate  a  variety  of  computer  systems  and  compare 
the  results.  NAVOCEANSYSCEN  supplied  ISA  with  their  Pre-Award  test  suite,  with 
instructions  to  also  include  certain  published  tests,  such  as  the  Dhrystone  and  Whetstone 
benchmarks.  It  was  ISA's  task  to  consolidate  all  the  designated  tests  into  one  easily- 
runnable  program  which  can  be  used  by  either  Government  or  Contractor  personnel  to 
test  all  computer  systems  that  are  candidates  for  procurement. 

NAVOCEANSYSCEN  selected  the  VAX  11/780  under  UNIX  4.3  BSD  as  the  baseline 
system  to  which  each  target  system  (procurement  candidates)  would  be  compared.  Our 
program,  bench,  collects  and  stores  test  results  from  the  baseline  system,  collects  and 
stores  results  from  all  the  target  systems,  and  produces  two  reports.  The  first  report 
compares  any  two  systems  which  the  user  selects,  and  the  second  report  summarizes  all 
the  test  data  in  one  report.  The  user  tells  the  program  which  of  these  tasks  he  wishes 
performed  through  the  use  of  options  on  the  input  line.  All  input  is  in  standard  UNIX 
format.  For  example, 

%  bench  -b 

would  execute  the  test  suite  and  store  the  results  as  the  baseline  data.  See  the  enclosed 
man  page  for  a  list  of  all  available  options. 

The  suite  includes  eighteen  tests  at  present.  The  user  is  able  to  specify  which  of 
the  eighteen  tests  will  be  used  by  modifying  an  ASCII  file  called  benchlist.  Benchlist 
includes  the  names  of  all  the  tests.  The  user  adds  or  deletes  the  comment  indicator  to 
tell  the  program  to  include  or  exclude  that  particular  test.  Additional  tests  may  be 
added  to  the  suite  by  modifying  the  benchlist  and  supplying  the  appropriate  code. 


Descriptive  Summary 

The  following  is  a  very  brief  description  of  the  tests  in  the  benchmark  suite. 


Fortran  Testa 


1.  Prime  Numbers 

This  program  generates  the  prime  numbers  from  0  to  8192  (optionally  printing  out 
the  results). 

2.  Calling  Sequence  and  Argument  Passing 

This  program  initializes  nine  variables,  passes  them  to  a  subroutine,  which  in  turn 
has  four  assignments.  This  sequence  is  repeated  one  million  times. 

3.  Random  Numbers 

This  program  tests  the  random  number  generator  by  generating  12,800  random 
numbers  and  checking  the  randomness. 

4.  Fast  Fourier  Transform 

This  program  performs  fast  fourier  transform  using  the  decimation-in-time  method 
(optionally  printing  out  the  results). 

5.  Matrix  Inversion 

This  program  performs  matrix  inversion  using  the  Gauss-Jordan  Reduction 
(optionally  printing  out  the  results). 

6.  Polynomial  Roots 

Roots  of  polynomials  are  calculated  using  the  Bairstow's  method  (optionally 
printing  out  the  results). 


Sieve  Tests 


7.  C  Sieve 

C  version  of  the  Sieve  of  Eratosthenes  prime  number  program. 

8.  Fortran  Sieve 

Fortran  version  of  the  Sieve  of  Eratosthenes  prime  number  program. 

9.  Pascal  Sieve 

Pascai  version  of  the  Sieve  of  Eratosthenes  prime  number  program. 


General  Tests 

10.  Whetstone 

A  C  version  of  the  original  Algol  benchmark,  "A  Synthetic  Benchmark"  by  H.  J. 
Cumow  and  B.  A.  Wichman.  Compiler  optimization  and  floating  point 
performance  are  tested. 

11.  Dhrystone 

This  program  contains  a  distribution  of  statements  which  are  considered  to  be 
representative:  53%  assignment,  32%  control  statements,  and  15%  procedure  and 
function  calls. 


Pre-Award  Testa 


12.  Block  Write 

This  program  creates  a  very  large  file  by  writing  8K  byte  blocks  one  after  the 
other. 

13.  Block  Read 

This  program  reads  the  file  created  by  the  block  write  program.  The  reads  are 
executed  in  8K  byte  blocks. 

14.  Sort 

A  shell  script  to  test  the  section  1  sort  call.  A  file  is  sorted  on  a  particular 
column  and  the  result  is  compared  to  a  presorted  file  to  test  the  results  of  the 
sort. 

15.  Integer  Arithmetic 

Addition,  subtraction,  multiplication,  and  division  are  performed  on  integer 
variables.  The  group  of  operations  is  executed  2.9  million  times. 

16.  Real  Arithmetic 

Similar  to  the  integer  arithmetic,  this  program  performs  addition,  subtraction, 
multiplication,  and  division  on  real  variables.  This  group  of  operations  is  executed 
600,000  times. 

17.  Large  Data  Space 

This  program  references  a  data  area  larger  than  real  memory  making  20,000 
references  to  random  locations. 

18.  Compile 

This  script  compiles  two  C  code  files  and  loads  the  two  object  files  into  a  single 
output  file. 
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Database  Format 

The  be.  chmark  test  data  is  stored  in  a  series  of  files  which  reside  in  the  current 
working  directory.  The  file  containing  the  baseline  data  is  called  baseline.  A  file  is 
created  for  each  system  tested  and  is  called  targetXXX,  where  XXX  is  a  3-digit  number 
assigned  by  the  benchmark  program  and  which  is  unique  to  each  system.  The  format  for 
both  the  baseline  file  and  all  the  target  files  is  identical. 

The  format  of  the  database  files  illustrated  on  the  following  page. 
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DEC  - 

VAX  1 1/780 

4  3  8SD  - 

5  •* - 

1  1  1  •+ - 

20  08  5.18  0  36 
20  04  5  09  0  25 
20  04  5  06  0  30 
0  51  0  05  0  10 

0  19  0.03  0  08 
0  50  0.05  0  07 
2  18  1.71  0  07 
1  99  1  67  0  06 
1  97  1  69  0  06 
1.71  1  33  0  05 
1  61  1  30  0  04 

1  61  1.30  0  07 

2  46  2  09  0  08 
2  36  2  08  0  05 
2  43  2  02  0  09 
10  05  5  26  0  38 
10  07  5  13  0  27 
10  09  5  24  0  47 

1 02  002  oo-; 

0  58  0  05  0  10 
1  28  0  01  0  14 
7  41  1  76  0  10 

5  20  1  69  0  04 
7  78  1  70  0  09 

6  25  1  35  0  08 

4  32  1  29  0  03 

6  56  1  37  0  08 

7  86  2  1 2  0  06 

5  73  2  10  0  05 

8  66  2  1  6  0.C7 
10  04  5  11  0  30 
10  08  5  09  0  28 
10  04  5  13  0  23 
0  71  0  02  0  09 
0  76  0  01  0  11 

0  76  0  07  0  06 

7  45  1  70  0  05 

6  49  1  69  0  03 
6  70  1  75  0  01 
5  49  1  35  0  03 

4  85  1  27  0  03 

5  03  1  30  0  04 

9  09  2  04  0  07 

8  40  2  1 2  0  06 
8  21  2  01  0  06 


EXECUTED 
ONE  AT  A 
TIME 


EXECUTED 
ALL  AT  ONCE 
(NO  LOAD) 


EXECUTED 
ALL  AT  ONCE 
WITH  EXTRA 
LOAD 


SYSTEM  MANUFACTURER 

SYSTEM  MODEL  NUMBER 

OPERATING  SYSTEM  VERSION 

NUMBER  OF  TESTS,  INCLUDING  ONE  FOR 

THE  COMPOSITE 

THREE  FLAGS  (  0  OFF,  1  ON)  FOR  THE  THREE 
WAYS  TO  EXECUTE  THE  SUITE:  SINGLY,  ALL 
AT  ONCE,  ALL  AT  ONCE  WITH  A  LOAD 

EACH  LINE  CONTAINS  THREE  NUMBERS: 

1.  REALTIME 

2.  USER  TIME 

3.  SYSTEM  TIME 

EACH  GROUP  OF  THREE  LINES  REPRESENTS 
THE  RESULTS  FROM  ONE  TEST  AND  EACH 
LINE  REPRESENTS  ONE  REPETITION  OF  THE 
TEST  THIS  IS  SET  IN  THE  C  CODE 
(CURRENTLY  REP  =  3) 


Reports 

Bench  produces  two  reports:  a  comparison  report  based  on  two  systems  of  the 
user's  choice,  and  a  summary  report  which  includes  all  systems  tested. 

The  comparison  report  is  invoked  when  the  user  specifies  the  -p  option  on  the 
command  line.  Bench  displays  a  list  of  those  systems  in  its  database,  and  prompts  the 
user  to  choose  two  systems  from  the  list.  The  comparison  report  displays  elapsed  time, 
user  time,  system  time,  and  percent  usage  for  each  test  and  each  system  chosen.  A 
composite  is  also  displayed.  The  composite  is  a  sum  of  all  systems  chosen  ar.d 
represented  as  if  it  were  a  separate  test.  The  elapsed  time  is  the  total  amount  of  time 
that  is  consumed;  the  "clock"  time.  The  user  time  is  the  amount  of  time  the  process 
spent  executing  nonprivileged  instructions  (e.g.,  arithmetic  calculations,  sorting, 
searching,  etc.).  System  time  is  the  time  the  process  spent  executing  privileged  (kernel) 
commands,  such  as  system  calls,  plus  some  system-level  overhead.  The  percent  usage  is 
that  portion  of  the  elapsed  time  that  is  actually  spent  executing  the  command.  It  is 
calculated  thusly: 


percent 

usage 


system  time  +  user  time 


elapsed  time 


x  100 


The  lowest  elasped  time  for  each  test  for  each  system  is  indicated  on  the  report 
by  an  asterisk  (*).  A  separate  column  is  displayed  for  the  elapsed  ratio.  The  first  figure 
in  the  elapsed  ratio  column  is  the  lowest  time  ratio,  which  is  the  ratio  of  the  lowest 
elapsed  time  of  the  second  system  to  the  lowest  elapsed  time  of  the  first  system  ,  or 


lowest  time  ratio 


lowest  elapsed  time  of  second  system 
lowest  elapsed  time  of  first  system 


The  second  figure  of  the  elasped  ratio  is  the  average  time  ratio.  The  average  time  ratio 
is  the  average  elapsed  time  of  the  second  system  divided  by  the  average  elapsed  time  of 
the  first  system,  or 


average  time  ratio  =  average  elapsed  time  of  second  system 

average  elapsed  time  of  first  system 

At  the  bottom  of  the  report,  an  average  elasped  ratio  is  given,  based  on  all  tests  except 
the  composite.  The  average  elapsed  ratio  is  the  average  of  each  lowest  time  ratio  (first 
figure)  and  the  average  of  each  average  time  ratio  (second  figure). 

The  summary  report  displays  elapsed  time  averages  for  all  systems  tested.  The 
test  results  are  normalized  to  the  baseline  system.  The  summary  report  is  invoked  when 
the  user  specifies  the  -e  option  on  the  command  line. 

Sample  output  for  both  comparison  report  and  the  summary  report  follow. 
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Ratios  are  displayed  as  the  [lowest  time  ratio,  average  time  ratio] 
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BRNCH(S) 


UNIX  PROGRAMMER'S  MANUAL 


BENCH(8) 


NAME 

bend)  -  benchmark  driver  and  result  comparison  generator 

SYNOPSIS 

bench  -b  [-▼]  [-a]  [-f]  [-1]  [-n  "make,  model,  version"] 

bench  -a  [-vj  f-a]  [-f]  [-1]  [-o  outfile]  (-n  "make,  model,  version"] 

bench  -p  [-a]  [-f]  [-1]  [-o  outfile ] 

bench  -e  [-a]  [-f]  [-1]  [-o  outfile ] 

DESCRIPTION 

Bench  is  a  benchmark  driver  program  to  time  the  execution  of  a  suite  of  tests  specified 
in  the  ASCII  file  benchlist.  The  defaults  to  bench  are  designed  to  allow  a  user  with 
little  or  no  understanding  of  the  options  to  establish  a  baseline  system  and  create 
comparisons  of  other  systems  with  the  baseline. 

The  user  can  override  the  defaults  by  using  the  options.  For  example:  data  can  be 
collected  without  generating  a  comparison  report;  an  output  filename  can  be  specified 
for  the  comparison  report;  a  summary  table  of  all  the  system  results  can  be  generated; 
and  selected  groups  of  tests  can  be  executed  without  running  the  complete  suite. 

The  available  flags  are: 

-b  Execute  the  test  suite  and  add  the  results  to  the  database  as  the  baseline  from 
which  comparisons  will  be  produced. 

-a  Execute  the  test  suite  and  add  the  results  to  the  database.  Unless  used  with  the  -o 
option,  no  comparison  report  will  be  generated. 

-p  Prepare  a  comparison  report  between  two  systems  of  the  user's  choice.  May  be 
used  with  the  -o  option;  default  is  to  the  line  printer. 

-e  Generate  a  summary  table  containing  normalized  elapsed  times  for  all  systems  in 
the  database.  If  no  system  has  been  assigned  as  the  baseline,  the  user  will  be 
prompted  for  a  system  to  use  as  a  baseline.  Default  is  to  the  line  printer. 

-v  Verbose:  causes  output  to  be  generated  to  standard  output.  This  information  is 
helpful  when  trying  to  follow  the  progress  of  the  driver.  Default  is  off. 

-n  "make,  model,  version" 

Use  the  make,  model,  and  version  of  the  system  to  identify  the  results.  This 
option  is  useful  when  executing  the  driver  in  a  batch  mode.  If  this  is  not  specified 
on  the  command  line,  the  user  will  be  prompted  for  make,  model,  and  version. 

'  -o  outfile 

Name  the  formatted  output  file  outfile.  By  default  the  output  file  is  created  by 
adding  the  last  three  digits  of  the  process  id  to  /tmp/bench. 


* 
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UNIX  PROGRAMMER'S  MANUAL 


BENCH(8) 


The  presence  of  any  of  the  -a,  -f,  and  -1  flags  cause  the  execution  to  be  limited  to  only 

what  is  specified.  (If  -a,  -f,  or  -1  are  not  specified,  the  default  sets  all  three  flags.) 

-a  This  flag  causes  tests  to  be  executed  one  after  the  other  with  no  extra  load  added 
to  the  system. 

-f  This  flag  causes  simultaneous  execution  of  the  tests  with  no  extra  load  added  to 
the  system. 

-1  This  flag  causes  simultaneous  execution  of  the  tests  with  extra  load  added  to  the 
system  at  the  same  time. 

EXAMPLES 

Execute  the  test  suite  on  the  current  machine  and  store  th<»  results  as  the  baseline. 


%  bench  -b 

Execute  the  test  suite  on  the  current  machine  and  store  the  results.  Also  input  the 
make,  model,  and  version  from  the  command  line. 

%  bench  -a  -n  "DEC,  VAX  8600 ,  4.3  BSD" 

Print  a  summary  of  current  database  into  outfile. 

%  bench  -e  -o  outfile 

FILES 


/t  mp/benchXXX 
baseline 
benchlist 
targetXXX 

BUGS 


formatted  output  of  the  comparison 

result  data  of  baseline  system 

the  path  of  test  and  the  printable  name 

result  data  of  system  to  be  compared  to  baseline 
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Instructions  for  Data  Collection  Using  the  Bench  Program 


Before  bench  can  be  used  on  any  system,  instructions  1  through  4  must  be  completed. 
All  tar  instructions  are  assuming  1600  bpi  on  drive  0. 

1.  Mount  tape  on  drive  0  at  1600  bpi. 

2.  Change  to  a  working  directory  with  at  least  3,000  blocks  free. 

3.  To  unload  the  tape,  type: 

%  tar  rv 

4.  To  compile  the  driver  program  and  test  suite,  type: 

%  make 

It  may  be  necessary  to  edit  the  bench.mk  file  to  alter  the  names  for  the  different 
compilers  with  optimizers  on. 

After  completing  steps  1-4  above,  any  of  the  remaining  sections  can  be  followed  to 
collect  data  or  display  previously  collected  data. 

To  establish  a  baseline  and  store  the  results  on  the  tape: 

1.  Type: 

%  bench  -b 

a.  Enter  the  make,  model,  and  version  of  the  system  when  prompted  by 
the  program. 

b.  Sit  back  and  relax. 

2.  To  store  the  baseline  file  on  the  tape,  type: 

%  tar  u  baseline 

3.  Remove  all  working  files  and  directories  from  the  disk  if  desired. 

To  add  a  target  system  to  the  database: 

1.  Type: 

%  bench  -a 

a.  Enter  the  system  description  when  prompted  by  the  program. 

b.  Sit  back  and  relax. 

2.  To  store  target  system  results  on  tape,  type: 

%  tar  u  target * 

If  a  print-out  is  desired,  skip  to  one  of  the  last  two  sections. 

3.  Remove  all  working  files  and  directories  if  desired. 


To  print  a  comparison  report  between  two  systems  in  the  database: 

1.  Type: 

%  bench  -p 

A  list  of  systems  in  the  database  will  appear  preceded  by  a  number.  The 
system  will  prompt  you  for  two  numbers  to  indicate  the  two  systems  to  be 
compared.  The  output  will  be  sent  to  the  line  printer. 

2.  Remove  all  working  files  and  directories  if  desired. 

To  print  a  summary  report  of  all  systems  in  the  database: 

1.  Type: 

%  bench  -e 

The  output  will  be  sent  to  the  line  printer. 

2.  Remove  all  working  files  and  directories  if  desired. 


lb 


Prim*  Numbers 


C 

C  prime. f 
C 

C  PROGRAM  TO  GENERATE  PRIME  NUMBERS 

C 

C  Compile  by:  fort  -O  prime. f  -o  prime 
C 

PROGRAM  PRIME 

C 

CONMOM/DAT/VALUE (8192) 

C 

C  INITILIZE  DATA  STRUCTURES 
C 

ILUM-6 

IPRT-0 

ICNT-512 

CUR-2. 

TOP- 3. 

1-1 

C 

C  CHECK  REMAINDER 

C 

5  I P ( AMOD ( TOP , CUR ) . EQ . 0 . ) GO  TO  10 

CUR-CUR* 1 . 

I F ( CUR . LT . TOP ) GO  TO  5 
C 

C  IP  WE  SCAN  PROM  2  THRU  TOP,  THEN  TOP  IS  A  PRIME  NUMBER 

C 

VALUE ( I ) -TOP 
1-1*1 
C 

C  SET  UP  POR  MEET  PRIME  NUMBER 

C 

10  TOP-TOP* 2. 

CUR-2. 

IP( I .LE. ICNTJGO  TO  5 
C 

C  PRINT  THE  PRIME  NUMBERS  WEVE  GENERATED 

C 

I P ( I PRT . EQ . 0 ) STOP 
DO  lb  1-1, ICWT.8 

WRITE ( ILUN , 9010 ) ( VALUE ( J ) ,J-I, 1*7) 

9010  FORMAT ( 8P10 . 0 ) 

15  CONTINUE 

STOP 

END 
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Calling  Sequence  and  Arguments  Passing 


calseq. f 

PROGRAM  TO  TEST  CALLING  SEQUENCE  AND  ARGUMENT  PASSING 
Compile  by:  fort  -O  calseq  -o  calseq 


PROGRAM  CALSEQ 
C 

Z-0. 

10  1-0 
J«1 

K«2 

L*3 

A-0. 

B-l. 

C-2. 

D-3. 

CALL  CALSEQKA,  I , B , J ,C , K , D, L ) 

Z-Z  +  l  . 

IF(Z.LT.l.E6)GO  TO  10 
STOP 

END 

C 

SUBROUTINE  CALSEQ1 ( A, I , B , J ,C , R, D, L ) 


I-J 

E-L 

RETURN 

END 


Random  Humber s 


C 

C  rndsk . f 
C 

C  PROGRAM  TO  PERFORM  A  CHECK  OF  THE  RANDOM  NUMBER 
C  GENERATOR  BT  PERFORMING  DIRECT  ACCESS  TO  A  DISK  FILE. 

C  THE  SUBROUTINE  WILL  USE  A  RANDOM  NUMBER  FROM  1  TO  256 
C  AS  THE  KEY  TO  READ  A  RECORD,  INCREMENT  THE  VALUE  READ, 

C  AND  WRITE  THE  NEW  VALUE. 

C 

C  Compile  by:  £ort  -O  rndi.sk. f  -o  rndsk 
C 

PROGRAM  RNDISK 

C 

ILUN-6 

IPRT-0 

ICNT-128 

FCNT-FLOAT< ICNT) 

FCHK-FCNT*100. 

IRAN-0 
isize  -  4 
B-rand( IRAN) 

OPEN ( ACCESS- ' DIRECT ' , 

1FILE- ‘ TEST ' , 

2 FORM- ' UNFORMATTED ' , 

C  3MAXREC-ICNT+1 , 

4RECL-isize*2, 

5STATUS- 'UNKNOWN ' , 

6 UNIT-4) 

C 

C  CREATE  FILE  WITH  EACH  RECORD  CONTAINING  ALL  ZEROS 

C 

DO  10  I -1 , ICNT 
IREC-I 

WRITE ( 4 , r ec - IREC ) FLOAT ( IREC) ,0. 

10  CONTINUE 

C 

C  GENERATE  ICNT* 100  RANDOM  NUMBERS 

C 

A-0 . 

20  IREC-IFIX ( FCMT* rand ( IRAN } ) >1 

IF(IREC.GE.l. AND. IREC. LE.ICNTJGO  TO  25 
WRITE ( ILUN, 9010 ) IREC 

9010  FORMAT ( '  RANDOM  NUMBER  OUT  OF  RANGE ’,16) 

25  I -IREC 

RSAD( 4 , r ec - I ) RNUM , COUNT 
COUNT -COUNT -t- 1 . 

I -IREC  ,, 

WRI TE ( 4 , r ec - I ) RNUM , COUNT 
A-A+l . 

I F ( A . LT . FCHK ) GO  TO  20 
C 

C  READ  PILE,  GET  MIN,  MAX  AND  AVERAGE  OF  RANDOM  NUMBER  GENERATOR 

C 

AMIN-9999. 

AMAX-0 . 


Random  Humber a 


30 


9000 


AVE-0 . 

DO  30  I-l.ICNT 
IREC-I 

RUAD( 4  ,  r  #C  - 1 REC ) RNUM ,  COUNT 

I F ( COUNT . GT . AMAX ) AMAX -COUNT 

IF ( COUNT , LT , AMIN ) AMIN -COUNT 

AVE-AVE+COUNT 

CONTINUE 

CLOSE( UNIT-4) 

AVE-AVE/FCNT 

IF( IPRT . EQ . 0 )STOP 

WRITE ( I LUN , 9000 ) AMIN, AMAX, AVE 

FORMAT (3FI 5.0) 

STOP 

END 
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Fast  Fourier  Transform 


fft.f 

PROGRAM  TO  PERFORM  A  FAST  FOURIER  TRANSFORM  USING  THE 
DECIMATION- IN-TIME  METHOD. 

Compile  by:  fort  -O  fft.f  -o  fft 

PROGRAM  FFT 

COMMOM/DAT/A( 4096 ) 

COMPLEX  A, U , W, T 

INITILIZE 

ILUM-6 
IPRT-0 

DO  25  LOOP-1, 10 
M-12 

ICNT-2**M 

PER-FLOAT( ICNT/16 ) 

PI-3.141592653589793 
DO  1  1-1 , ICNT 
B-SIH( 2. *PI*FLOAT( I)/PER) 

A ( I ) -CMPLX ( B , 0  .  ) 

1  CONTINUE 

N«2**M 
NV2-N/2 
NMl-N-1 
J-l 

DO  7  I  -1 ,  NM1 

IF( I . GE. J)GO  TO  5 

T-A(J) 

A(  J ) »A( I  ) 

A  (  I ) -T 

5  K-NV2 

6  IF( K . GE . J )GO  TO  7 

J-J-K 

K-K/2 
GO  TO  6 

7  J-J+K 
PI-3.141592653589793 
DO  20  L-1,M 
LE-2**L 

LE1-LE/2 
U-(l. ,0. ) 

W-CMPLX ( COS ( P I /FLOAT ( LEI ) ) , SIN( PI /FLOAT ( LEI ) ) ) 
DO  20  J-l, LEI 
DO  10  I»J,N,LE 
IP-I+LE1 
T»A( IP)*U 
A(  IP) -A( I ) -T 
10  A(I)-A(I)+T 

20  U-U*N 

25  CONTINUE 

IF( I PRT . EQ . 0 ) STOP 
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Fast  Fourier  Transform 


DO  30  I«l, 128,4 

WRITE(  ILUN,  9000)  (A(J)  ,J*I,  1  +  3) 
9000  FORMAT (4G1 5. 6) 

30  CONTINUE 

STOP 
END 
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Matrix  Inversion 


matrix. £ 

MATRIX  INVERSION  USING  GAUSS -JORDAN  REDUCTION 
INVERTED  MATRIX  OVERLAYS  ORIGINAL  MATRIX  IN  MEMORY 
PARITAL  PIVOTING  IS  NOT  USED 

Compile  by;  fort  -O  matrix  -o  matrix 

PROGRAM  MATRIX 

COMMON/DAT/A (15,15) 

DOUBLE  PRECISION  A 

ILUM-6 

IPRT-0 

DO  10  LOOP-1,10000 

N»4 

A( 1 , 1 ) “1 . 

A<1,2)“1. 

A( 1 , 3 ) “1 . 

A(l,4)-1. 

A(  2, 1 ) *4 . 

A(  2 , 2 ) *5 . 

A(  2 , 3 ) *6  . 

A{ 2 , 4 ) *7 . 

A<3,1)“6. 

A(  3 , 2 ) “10 . 

A( 3 , 3 ) *15 . 

A( 3 , 4 ) *21 . 

A( 4 , 1 ) “12 . 

A( 4 , 2 ) *30 . 

A(  4 , 3 ) “60 . 

A( 4 , 4 ) “105 . 

CALCULATE  ELEMENTS  OF  REDUCED  MATRIX 
DO  6  K-1,R 

CALCULATE  NEW  ELEMENTS  OF  PIVOT  ROW 

DO  4  J-l ,H 

IF  (  J  .  EQ  .  K  )GO  TO  4 

A(K, J)*A(K, J)/A(K,K) 

CONTINUE 

CALCULATE  ELEMENT  REPLACING  PIVOT  ELEMENT 
A(K,K)“1./A(K,K) 

CALCULATE  NEW  ELEMENTS  NOT  IN  PIVOT  ROW  OR  PIVOT  COLUMN 

DO  5  I -1 , N 
I F  ( I .  EQ .  K  )  GO  TO  5 
DO  5  J-l ,N 
I F ( J . EQ . K ) GO  TO  5 
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Matrix  Inversion 


A(I,J)*A(I,J)-A(K,J)*A(I,K) 

CONTINUE 

CALUCLATE  REPLACEMENT  ELEMENTS  FOR  PIVOT  COLUMN-EXCEPT  PIVOT  ELEMENT 

DO  6  1*1, N 

IF ( I . EQ, K ) GO  TO  6 

A( I,K)*-A( I , K ) *A( K , K ) 

CONTINUE 
I  CONTINUE 

OUTPUT  INVERTED  MATRIX 

IF( I PRT . EQ . 0 ) STOP 

WRITE ( ILUN, 8 ) ( (A(I,J),J=1,N),I*1,N) 

FORMAT ( 4F1 6 . 4 ) 

STOP 

END 
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Polynomial  Roots 


C 

C  roots. f 
C 

C  ROOTS  OP  POLYNOMIAL  8Y  BAIRSTOWS  METHOD 
C 

C  Compile  by:  fort  -O  roots. f  -o  roots 
C 

PROGRAM  ROOTS 
C 

DIMENSION  A(30) ,8(30) ,C(30) 

ILUN-6 

IPRT-0 

IP( IPRT.EQ.O)GO  TO  200 

IPRT-0 

JPRT-1 

200  DO  100  LOOP-1 , 10000 

If ( LOOP. HE. 10000 )GO  TO  220 
I P ( JPRT . EQ . 0 ) GO  TO  220 

IPRT-1 

220  UI-0. 

VI-0. 

EPS I - 1 . E - 6 
N-5 

A(l)— 3. 

A(2)— 10. 

A( 3 ) >10 . 

A(  4 ) -44 . 

A( 5 ) -48 . 

C 

C  SEE  IP  H-0,1,  OR  GREATER  THAN  1 
C 

40  IP(N-1 ) 100 ,5,7 

5  P— A(l) 

Q-0 . 

IT-1 

IP( IPRT .NE. 0 ) WRITE ( ILUN , 6 )N , P, Q , IT 

6  FORMAT ( '  X(',I2,')  - ' , 2X , P8 . 4 , 6X , P8 . 4 , 10X , I 3 ) 
GO  TO  100 

C 

C  SEE  IP  N-2  OR  IP  N  IS  GREATER  THEM  2 

C 

7  I P  ( M .  EQ .  2 )  GO  TO  8 

GO  TO  13 

8  U-A(l) 

V-A( 2 ) 

IT-1 

9  P--U/2. 

RAD-U**  2-4 . *V 

C 

C  CHECK  THE  SIGN  OP  U**2-4.*V 
C 

IP( RAD.GT . 0 . )GO  TO  12 
RAD- -RAD 
Q-SQRT(RAD)/2. 

IP( IPRT. NE. 0) WRITE ( ILUN , 6 ) N , P, Q , IT 


Polynomial  Roots 


M-M-l 

Q.-Q 

90  IF{  IPRT. ME.  0 } WRITE ( ILUN ,  6 ) M, P,Q , IT 

10  M-M-l 
C 

C  CHECK  TO  SEE  IP  N  IS  GREATER  THEM  ZERO 
C 

IP ( M . LE . 0 ) GO  TO  100 
DO  11  I*1,M 

11  A(I)«B(I) 

GO  TO  40 

12  Q-SQRT(RAD)/2. 

W-P 

Z-Q 
P-P+Q 
Q-0 . 

IP( IPRT. ME. 0 ) WRITE ( ILUN, 6)M,P,Q, IT 

M-M-l 

P-W-Z 

GO  TO  90 

13  U-UI 
V-VI 
IT-1 

C 

C  CALCULATE  THE  B  VALUES 
C 

50  B  ( 1 ) >A( 1 ) -U 

B( 2) «A( 2 ) -B( 1 ) *U-V 
DO  14  K-3  ,N 

14  B(K)-A(K)-B(E-1)*U-B(K-2)*V 
C 

C  CALCULATE  THE  C  VALUES 
C 

C(l)-B(l)-U 

C(2)-B(2)-C(1)*U-V 

M-M-l 

DO  15  K-3 ,M 

15  C(K)-B(K)-C(K-1)*U-C(K-2)*V 
C 

C  CALCULATE  DELU  AMD  DELV 
C 

IF(M.GT. 3 )GO  TO  17 
DEMOM-C (M-l)-C(M-2)**2 
I P ( DEMON . EQ . 0 . ) GO  TO  30 
DELU-( B(M ) -B(M-l)*C(M-2) ) /DEMON 

16  DELV-(C(H~1 )*B(M-1)-C{M-2)*B(M) ) /DEMOM 
GO  TO  18 

17  DEMOM-C (M-l )*C(M-3)-C(M-2)**2 
IP ( DEMOM. EQ . 0 )GO  TO  30 

DELU-(B(M ) *C(M-3)-B(M-l)*C(M-2) ) /DEMOM 
GO  TO  16 
C 

C  CALCULATE  MEW  U  AMD  V  VALUES 
C 
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Polynomial  Roots 


18  U-U+DELU 
V-V+DELV 

SUN -AES ( DELU ) +ABS ( DELV ) 

C 

C  STORE  THE  FIRST  SUM  CALCULATED 
C 

IF( IT. EQ. 1 )GO  TO  19 
GO  TO  20 

19  STORE-SUN 
GO  TO  21 

20  IF( IT. EQ . 50)GO  TO  28 
IF( IT.GE. 200) GO  TO  26 

21  I P ( SUM . LE . EPS  I ) GO  TO  9 
IF( IT . EQ . 100 )GO  TO  23 

22  IT-IT+1 
GO  TO  50 

23  IP( IPRT.NE. 0 )WRITE( ILUN , 24 ) 

24  FORMAT ( '  CONVERGENCE  IS  SLOW') 

IF ( IPRT.HE. 0) WRITE ( ILUN,25)U,V 

25  FORMAT (  '  U-',E14.7,'  V-',E14.7) 

GO  TO  22 

26  IF( IPRT . ME . 0 ) WRITE ( ILUN ,  27 ) 

27  FORMAT ( '  STOPPED  AFTER  200  ITERATIONS') 

IF( IPRT . NE . 0 ) WRITE ( ILUN,25)U,V 

GO  TO  100 
C 

C  SEE  IF  SUM  AFTER  50  ITERATIONS  EXCEEDS  FIRST  SUM  STORED 

C 

28  I F ( SUM . LT . STORE ) GO  TO  21 
IF( IPRT. NE.0)WRITE( ILUN ,29) 

29  FORMAT ( '  DIVERGENCE  OCCURRING') 

IF< IPRT. ME.0)WRITE< ILUN, 25 )U,V 
GO  TO  100 

30  IF( IPRT.NE. 0)WRITE( ILUN, 31) 

31  FORMAT ( '  DENOMINATOR  IS  ZERO') 

GO  TO  100 

100  CONTINUE 

STOP 
END 


C  Sieve 


/* 

•  sieve. c 

• 

*  Eratosthenes  Sieve  Prime  Humber  Program  in  C  */ 

* 

*  Compile  by;  cc  -O  sieve. c  -o  csieve 

* 

v 

♦define  true  1 
♦define  false  0 
♦define  size  8190 

char  flags(size  +  1]; 


main( )  { 

int  i .prime, k, count, iter ; 

printf  (  "100  iterationsm"  ) ; 
for(iter  ■  l;iter  <*  100;iter  ++)  { 
count*0 ; 

for(i  *  Q;i  <«  size;i  ++ ) 
flags[i]  *  true; 
for(i  *  0;i  <*  size;i  ♦+)  { 
i f ( f lags  { i J  )  { 

prime  ■  i  ♦  i  ♦  3; 
for ( k-i+pr iae ; k<-s ize ; k+«pr ime ) 
flags [k]  «  false; 

count++; 


printf ("%d  is  largest  of  %d  pr imes . »n" , pr ime , count ) ; 


Fortran  Sieve 


c 

c  sieve. f 
c 

c  eratosthenes  sieve  with  Rnuth's  optimization 
c 

c  Compile  by:  fort  -O  sieve. f  -o  fsieve 
c 

integer  i , j , k , iter , pr ime, count 
logical  flags( 8191 ), last 

wr i te( 6 , 10 ) 

10  format  ('  100  iterations') 
do  20  iter  *  1,  100 

count  *  0 
do  30  i  -  1,  8191 
30  flags(  i  )  *  . true, 

last  *  .false, 
do  40  i  *  1,  8191 

if  (.not.  flags(i))  go  to  50 
prime  *  i  +  i  +  1 
count  ■  count  +  1 
c  write(6,ll)  prime 

11  format  (lx,i6) 

if  (last)  go  to  50 

k  *  ( pr ime*pr ime  -  1)  /  2 
c  k  *  i  +  prime 

do  60  j  -  k,  8191,  prime 
60  flags(j)  ■  .false. 

if  (prime  . ge .  127)  last  ■ 
50  continue 

40  continue 

20  continue 

wr i te ( 6 , 12 )  count 

12  format  (lx,  i6,  '  primes') 
end 


. true . 


Pascal  Sieve 


( *  sieve . p  * ) 

(•  Eratosthenes  Sieve  Prime  Humber  Program  in  Pascal  *) 
(•  Compile  by.  pi  sieve. p  *) 
program  pr ime( output ) ; 
cons  t 

size  *  8190; 
var 

flags  :  array  [0..size]  of  boolean; 

1 , pr ime , k , cnt , i ter  :  integer; 

begin 

writeln('100  iterations’); 
for  iter  :»  1  to  100  do  begin 
cnt  : *  0 ; 

for  i  : =  0  to  size  do 

flags [i]  true; 
for  i  : *  0  to  size  do 

if  flags (i]  then  begin 
prime  :=  i+i+3; 
k  :*  i  ♦  prime; 
while  k  <*  size  do  begin 

flags[k]  : -  false; 
k  :■  k  +  prime 
end  ; 

cnt  :■  cnt  *  1 
end; 

end; 

wr i teln( cnt ,  '  primes') 
end. 


K) 


Whet  stone 


whet  .  c 

Whetstone  benchmark  in  C.  This  program  is  a  translation  of  the 
original  Algol  version  in  " A  Synthetic  Benchmark*  by  H.J.  Curnow 
and  B.A.  Wichman  in  Computer  Journal,  Vol  19  |1,  February  1976. 

Used  to  test  compiler  optimisation  and  floating  point  performance. 

Compile  by:  cc  -0  -s  -o  whet  whet.c 

or:  cc  -O  -DPOUT  -s  -o  whet  whet  c 

if  output  is  desired. 


♦define  iTERATIOWS  10  /*  1  Million  Whetstone  instructions  */ 


♦include  "math.h* 


double 

double 

i  r  t 


Ml,  m2,  m3,  m4,  m,  y,  z,  t,  tl,  1 2 ; 
el[4]; 

l,  'j  ,  k,  1,  nl,  n2,  n3,  n4,  n6,  n7,  n8,  n9,  nlC,  nil; 


main( j 


/*  initialize  constants  */ 

t  -  0.499975; 

tl  •  0.50025; 

1 2  -  2.0; 

/•  set  values  of  module  weights  */ 


nl 

tt 

0 

* 

ITERATIOWS 

n2 

W 

12 

• 

ITERATIOWS 

n3 

m 

14 

* 

ITERATIOWS 

n4 

m 

345 

• 

ITERATIOWS 

n6 

m 

210 

* 

ITERATIOWS 

n  7 

m 

32 

• 

ITERATIOWS 

n8 

m 

899 

• 

ITERATIOWS 

n9 

m 

616 

• 

ITERATIOWS 

nlO 

m 

0 

* 

ITERATIOWS 

nil 

m 

93 

• 

ITERATIOWS 

/*  MODULE  1:  simple  identifiers  •/ 
Ml-  1.0; 

m2  •  m3  -  m 4  »  -1.0; 


1; 

1 

»  n  1  ; 

1 

■f  » 

1  )  ' 

Ml 

* 

Ml  ♦ 

m2 

•f 

x  3  m  4 

* 

t 

m2 

* 

Ml  ♦ 

m2 

- 

m3  m4 

* 

t 

m3 

* 

Ml  - 

m2 

♦ 

x  3  ♦  x  4 

* 

t 

m4 

* 

-Ml  + 

m2 

♦ 

m3  ♦  M  4 

* 

t 

f  or  (  i 


Whetstone 


I 1 fdef  POUT 

pout(nl,  nl,  nl,  si,  «2,  ■ 3 ,  «4  )  ; 

fend  1  f 


/*  MODULE  2:  array  element*  •/ 


el  [ 0 1 

«  1.0, 

el  [  1 1 

-  el [2] 

*  el  [  3  1  • 

-1.0; 

for  ( 

l  «  1;  l 

<  «  n  2  ;  l 

4*1)  ( 

el  10] 

-  (  el(0] 

♦  el  11] 

4  el ( 2]  -  el [ 3] 

)  *  t 

el[l! 

*  (  el (0] 

♦  e  1  1 1  ] 

-  el [ 2 ]  4  el [ 3] 

)  *  t 

el  (2) 

-  <  el  10] 

-  el [1] 

4  el[ 2]  4  el ( 3 ] 

)  *  t 

el  [  3] 

*  ( -el  [0] 

♦  el  ]  1  ] 

4  el [ 2]  4  el [ 3 ] 

)  *  t 

} 

I i fdef  POUT 

po  u  t  ( n  2 ,  n3 ,  n2,  el[0],  el[l),  e 1  [  2 ]  ,  e 1 [ 3 ]  ) ; 

tend  if 

/*  MODULE  3:  array  a*  parameter  •/ 

for  (i  »  1;  i  < ■  n3;  i  4«  1) 
pa( el ) ; 

•  i fdef  POUT 

poat(n3,  n2 ,  n2,  el[0],  el[l],  el(2],  e 1 [ 3 ]  )  ; 

tend  l  f 

/•  MODULE  4:  conditional  jump*  */ 


J  '  , 

for  (i  *  1;  i  <■  n4;  i  ♦*  1)  { 

if  (j  --  1) 

j  -  2; 

else 

j-3; 

if  ( j  >  2) 

j  -  0; 

else 

3-i; 

if  ( j  <  1  ) 

j  -  i; 

else 

3  -  0; 

) 

tifdef  POUT 

pout ( n4 ,  ],  3,  *1,  s  2 ,  »3,  *4); 


tend  i  f 

/•  MODULE 

5; 

omitted  */ 

/•  MODULE 

fl  : 

integer  arithmetic  •/ 

Whetstone 


j  -  l: 

w  -  2 ; 

1  -  3; 

for  (i  »  1;  i  < ■  n6;  i  +  *  1)  { 

j-j*(k-j)*(l-k>; 
k  -  1  *  k  -  (1  -  j)  *  k; 

1  -  (1  -  k)  *  (k  ♦  j); 

el [ 1  -  2)  «  j  ♦  k  +  1 ;  /*  C  arrays  are  zero  based  */ 

e  1  [  k  -  2  ]  *  j  *  k  •  1 ; 

#ifdef  POUT 

po  u  t ( n  6  ,  j,  k,  el  [  0 ]  ,  el(l],  el[2],  el[3]); 

•  end  i  f 

/•  MODULE  7:  trig,  functions  */ 
s  *  y  *  0. 5; 

for(i  *  1;  i  <*  n7;  i  +*1)  { 

x  ■  t  *  atan( t2*s in( x ) *cos ( x )/( cos ( x+y ) +cos( x-y ) -1 . 0 ) ) ; 
y  -  t  *  atan( 1 2*s in( y ) *cos ( y )/( cos ( x+y ) +cos( x-y ) -1 . 0 )  )  ; 

} 

lifdef  POUT 

pout ( n7 ,  j,  k,  x,  x,  y,  y); 

lend! f 

/•  MODULE  8:  procedure  calls  */ 

x»y«*«l.Q; 

for  (i  •  1;  i  <«  n8;  i  +»1) 
p3 ( x ,  y,  fcz  )  ; 

tifdef  POUT 

pout ( n8 ,  j,  k,  x,  y,  z,  z); 

lend  i  f 

/*  MODULE 9 :  array  references  */ 

j  ■  1; 

k  *  2; 

1  *  3; 


e l ( 0 ]  -  10; 
el  111  -  2.0; 
el [2]  -  3.0; 

for(i  *  1;  i  <-  n9;  i  +•  1) 

p0  (  )  ; 

I i fdef  POUT 

pout ( n9 ,  3,  k,  el [0] ,  e 1 [ 1 ] ,  el[2],  el[3]); 


lend  i  f 


Whetstone 


/*  MODULEIO:  integer  arithmetic  */ 

j  *  2; 

k  «  3; 


for ( i  ■ 

i; 

i  <* 

nlO; 

i  -t-1) 

{ 

j 

*  j  ♦ 

k; 

k 

-  j  + 

k; 

j 

*  k  - 

j  ; 

} 

k 

»  k  - 

j  - 

j; 

POUT 

pout( nlO 

9 

j  *  ^  9 

*i. 

x2,  x3. 

x4  )  ; 

•endif 

/*  WODULE11:  standard  functions  */ 
x  -  0.75; 

for(i  «  1;  i  <■  nil;  i  +*1) 

*  -  sqr t (  exp(  log(x)  /  tl)); 


• i fdef  POUT 

pout( nil ,  j,  k,  x,  x,  x,  x); 

♦endif 
exit  (0); 

} 


pa(e) 

double  e[4]  ; 

register  int  j; 
j  -  0; 

lab : 

•  10)  -  (  e [0]  +  ell]  +  e 1 2]  -  e[3]  )  *  t; 

*U]  -  (  e  10]  +  ell]  -  e[2]  +  e[3]  )  *  t; 

e l 2]  -  (  *10]  -  e ( 1]  +  e l 2]  +  *13]  )  •  t; 

*[3]  «  (  -e[0]  ♦  ell]  +  e l 2]  ♦  e(3J  )  /  t2; 

j  i; 
if  (j  <  6) 

goto  lab; 


p3(x,  y,  z) 
double  x,  y, 

x  * 

y  - 

*Z  * 

} 


*s; 

t  *  ( *  +  y ) ; 
t  *  (x  f  y); 
(x  ♦  y)  /t2; 


U 


Whetstone 


ell j]  -  ellk] ; 
el [k]  -  el [ 1] ; 
•  Ml]  ■  •lljl; 


|i  fdef  POUT 

pout ( n,  j,  k,  xl,  x2,  x3 ,  x4 ) 

int  n,  j,  k; 

double  xl,  x2,  x3,  x4; 

{ 

printf("%6d%6d%6d  %5e  «5e  «5e  »5em" 

n,  j,  k,  xl,  x2,  x3 ,  x4); 

lenda  f 


Dhrystone 


dry.c 


"DHRYSTONE"  Benchmark  Program 


Version:  C/1.1,  12/01/84 


Date : 
Author : 


Compile : 


Run : 


PROGRAM  updated  01/06/86,  RESULTS  updated  02/17/86 

Reinhold  P.  Weicker,  CACM  Vol  27,  No  10,  10/84  pg.  101 
Translated  from  ADA  by  Rick  Richardson 
Every  method  to  preserve  ADA-likeness  has  been  used, 
at  the  expense  of  C-ness. 


cc  -O  dry.c  -o  drynr  :  No  registers 

cc  -O  -DREG=register  dry.c  -o  dryr  :  Registers 

drynr;  dryr 


The  following  program  contains  statements  of  a  high-level  programming 
language  (C)  in  a  distribution  considered  representative: 


assignments  53% 
control  statements  32% 
procedure,  function  calls  15% 


100  statements  are  dynamically  executed.  The  program  is  balanced  with 
respect  to  the  three  aspects: 

-  statement  type 

-  operand  type  (for  simple  data  types) 

-  operand  access 

operand  global,  local,  parameter,  or  constant. 


The  combination  of  these  three  aspects  is  balanced  only  approximately. 


The  program  does  not  compute  anything  meaningfuJ.1,  but  it  is 
syntactically  and  semantically  correct. 


U- 


Dhrystone 


/*  Accuracy  of  timings  and  human  fatigue  controlled  by  next  two  lines  */ 
/•♦define  LOOPS  50000  /*  Use  this  for  slow  or  16  bit  machines  */ 

tdefine  LOOPS  500000  /•  Use  this  for  faster  machines  */ 

/*  Compiler  dependent  options  */ 

iundef  NOENUM  /*  Define  if  compiler  has  no  enum's  */ 

iundef  NOSTRUCTASSIGN  /*  Define  if  compiler  can't  assign  structures  */ 

/*  define  only  one  of  the  next  two  defines  */ 

♦define  TIMES  /*  Use  times(2)  time  function  */ 

/•♦define  TIME  /*  Use  time(2)  time  function  */ 

/*  define  the  granularity  of  your  times(2)  function  (when  used)  */ 

♦define  HZ  60  /*  times(2)  returns  1/60  second  (most)  */ 

/•♦define  HZ  100  /*  times(2)  returns  1/100  second  (WECo)  */ 

/*  for  compatibility  with  goofed  up  version  */ 

/•♦define  GOOF  /*  Define  if  you  want  the  goofed  up  version  */ 

♦ ifdef  GOOF 

char  VersionI]  *  "1.0"; 

♦  else 

char  Version[]  *  "1.1"; 

♦endi f 

♦ifdef  NOSTRUCTASSIGN 

♦define  structassignld,  s)  memcpy(6(d),  Ms),  sizeof(d)) 

♦  else 

♦define  structassign(d,  s)  d  *  s 

♦  end  if 

♦ifdef  NOENUM 

♦define  Identl  1 

♦define  Ident2  2 

♦define  Ident3  3 

♦define  Ident4  4 

♦define  Ident5  5 

typedef  int  Enumeration; 

♦  else 

typedef  enum  {identl,  Ident2,  Ident3,  Ident4,  Ident5}  Enumeration; 

♦endi f 

typedef  int  OneToThirty; 
typedef  int  OneToFifty; 

typedef  char  Capi talLetter ; 
typedef  char  Str ing30 [ 31 )  ; 
typedef  int  ArraylDim[51] ; 

typedef  int  Ar ray2Dim [51]  [51]  ; 

struct  Record 

{ 


struct  Record 

Enumeration 

Enumeration 


•PtrComp; 
D  i  s  c  r  ; 
EnumComp ; 


Dhrystone 


OneToFifty  IntComp; 

String 30  StringComp; 


typedef  struct  Record  RecordType; 

typedef  RecordType  *  RecordPtr; 

typedef  int  boolean; 

•define  NULL  0 

•define  TRUE  1 

•define  FALSE  0 

•ifndef  REG 
•define  REG 
•endif 

extern  Enumeration  Funcl{); 

extern  boolean  Func2(); 

tifdef  TINES 
•include  <sys/types .h> 

•include  <sys/times .h> 

•endif 

main( ) 

ProcO( ) ; 
exit{0) ; 


/* 

*  Package  1 
*/ 

int 

boolean 

char 

char 

Ar raylOim 
Array 2Dim 
RecordPtr 
RecordPtr 


IntGlob; 
BoolGlob; 
Char IGlob; 
Char 2Glob; 
Array IGlob; 
Array2Glob; 
PtrGlb; 
PtrGlbNext ; 


Proc0( ) 

{ 

OneToFifty 
REG  OneToFifty 
OneToFifty 
REG  char 
REG  char 
Enumeration 
Str ing30 
Str ing30 
extern  char 


ZntLocl ; 
IntLoc2; 
IntLoc3; 
CharLoc; 
Char Index; 
EnumLoc ; 

Str inglLoc ; 
Str ing2Loc; 
•malloc { ) ; 


Dhr  y s  tone 


tifdef  TIME 
long 
long 
long 
long 

register  unsigned  int 


time(  ) ; 
star ttime ; 
bencht ime; 
null time; 
i  ; 


starttime  *  time(  (long  *)  0); 
for  (i  *  0;  i  <  LOOPS;  +  +i); 

nulltime  *  time(  (long  *)  0)  -  starttime;  /*  Computes  o ’ head  of  loop  */ 

fendif 

tifdef  TIMES 

t ime_t 
t ime_t 
t ime_t 
struct  tms 

register  unsigned  int 


s  tar  1 1 ime ; 
bencht ime ; 
nulltime; 
tms  ; 
i ; 


times(itms);  starttime  *  tms . tms_ut ime ; 
for  ( i  *  0 ;  i<  LOOPS;  ++i); 
t imes ( 6 tms ) ; 

nulltime  *  tms . tms_ut ime  -  starttime;  /*  Computes  overhead  of  looping  */ 

fendif 


PtrGlbHext  *  (RecordPtr)  mal loc( si zeof ( RecordType ) ) ; 

PtrGlb  »  (RecordPtr)  malloc ( s izeof ( RecordType )) ; 

PtrGlb->PtrComp  -  PtrGlbHext; 

PtrGlb->Discr  *  Identl; 

PtrGlb->EnumComp  *  Ident3; 

PtrGlb->IntComp  *  40; 

strcpy( PtrGlb- >Str ingComp,  " DHRYSTONE  PROGRAM,  SOME  STRING* ) ; 
tifndef  GOOF 

s tr cpy ( Str inglLoc ,  "DHRYSTONE  PROGRAM,  l’ST  STRING");  /"GOOF*/ 

tend  i  f 

Ar ray 2Glob [ 8 ]  [7]  »  10;  /*  Was  missing  in  published  program  */ 


/**♦***“*• . . 

--  Start  Timer  -- 

. . 

tifdef  TIME 

star ttime  «  time(  (long  *)  0); 

fendif 

tifdef  TIMES 

t imes ( 6 tms ) ;  starttime  -  tms. tms  utime; 

fendif 

for  (i  ■  0;  i  <  LOOPS;  ++i) 


Proc5( ) ; 

Proc4 (  )  ; 

IntLocl  »  2; 

IntLoc2  »  3; 

Strcpy (String 2 Loc,  "DHRYSTONE  PROGRAM,  2'ND  STRING"); 
EnumLoc  «  Ident2; 


V) 


Dhrystone 


BoolGlob  *  !  Func2(Str inglLoc,  String2Loc); 
while  (IntLocl  <  IntLoc2) 

{ 

IntLocl  *  5  *  IntLocl  -  IntLoc2; 

Proc7  (  IntLocl ,  l'ntLoc2,  fclntLod): 

++IntLocl : 

) 

Proc8(ArraylGlob,  Array2Giob,  IntLocl,  IntLoc3); 

Prod  ( PtrGlb )  ; 

for  (Charlndex  *  ’A1;  Charlndex  < »  Char2Glob;  ++Charlndex) 
if  (EnumLoc  **  Funci(CharIndex,  *C')) 

Proc6 ( Identl ,  fcEnumLoc); 

IntLoc3  *  IntLoc2  *  IntLocl; 

IntLoc2  *  IntLocl  /  IntLocl; 

IntLoc2  ■  7  *  (IntLoc3  -  IntLoc2)  -  IntLocl; 

Proc2( tint Loci )  ; 


/***““ . ***** 

--  Stop  Timer  -- 

* . 


t i fdef 


lend  i  f 
I l fdef 


lendif 


TIME 

benchtime  «  time(  (long  *)  0)  -  start  time  -  null  time; 
pr  intf  (  "Dhry»tone(  %s )  time  for  ftld  passes  «  Udm", 

Version, 

(long)  LOOPS,  benchtime); 

pr intf ( "This  machine  benchmarks  at  %ld  dhrystones/second'n" , 
((long)  LOOPS)  /  ►'enchtime); 

TIMES 

times(fctms) ; 

benchtime  »  tms .  tn>i_utime  -  starttime  -  nulltime; 
pr intf ( "Dhrystone ( %s)  time  for  %ld  passes  «  %ld*n", 

Version, 

(long)  LOOPS,  bencht ime/HZ ) ; 

pr intf ("This  machine  benchmarks  at  %ld  dhrystones/second*n" , 
((long)  LOOPS)  *  HZ  /  benchtime); 


) 


Prod  (  PtrPar  In ) 

REG  RecordPtr  PtrParln; 

{ 

♦define  NextRecord  (*( PtrPar In- >PtrComp) ) 

s t rue tass i gn ( Hex t Record ,  "PtrGlb); 
PtrParIn->IntComp  «  5; 

NextRecord . IntComp  -  PtrParIn->IntComp; 
NextRecord . PtrComp  «  PtrPar In->PtrComp; 
Proc3(NextRecord  PtrComp); 
if  (  Nex  t  Peco  rd  .  Di  sc  r  **  Ident.J) 

{ 


40 


Oh  f  yilon* 


Kt«t  R»cord .  IntComp  -  h , 

Pr  or6(  Pt  r  P»i  In  'EnumComp,  lN*a »  Rucoi  cl  .  EnnmComp  )  ; 
Ea«t Ruroi d . Pt rComp  •  PtrGlb  -Pt rComp; 

^  Prorf(R»«tR*cnrd  IntComp,  10,  tR«a t Rfcn i d . I n t Comp ) 

•  1 

it  r  or  t  a  •  ■  i  q  n  (  *  P  t  rPti  In,  Rail  Rirnt  rl )  ; 
fundif  R*«tR*rmd 

I 


ProrJj IntPir IO) 

OnaToF  iffy  *  I  n  t  Pa  r  If); 

( 

REG  OnaToFifty  Inti. or; 

REG  Enumarat  ion  EnumLoc  ; 


I 


I  n  t  Lot" 
for  (  :  ; 
( 


I 


-  * Int  Par  IO  ♦  10; 


i f  (Char IGlob  -»  '  A  1  ) 

( 

I  n  t  l.or  ; 

MntParlO  -  IntLor 
EnumLoc  »  Idant  1  ; 

} 

if  (EnumLoc  Idant  I  ) 
hraak  ; 


I  n  t  G 1  o  h ; 


Pr  oc I ( Pt  r ParOut  ) 

RarordPt  r  ‘Pt  rPirOut  ; 

( 

if  ( Pt  rGl h  i -  HULL ) 

•PtrParOut  -  PtrGlb  -Pt  rComp; 

a  1  sa 

IntGlob  -  100; 

Pr  oc  1(10,  IntGlob,  API rC 1 b  -IntComp); 


Pr  or  4  (  ) 

( 

REG  boolaan  Rooll.oc; 

Ron  1  Lor  -  Char IGlob  »•  ’A’; 
Rooll.oc  |»  Roo  IGlob; 

Char  JIG  I  oh  -  '  R  '  ; 

I 


Pr  oc'.  (  ) 

I 

Cha  r  I  G  1  of*  -  '  A  '  ; 

Roo I G I  oh  -  FALSE ; 


1  I 


Dhrystone 


} 

extern  boolean  Func3(); 

Proc6(EnumParIn,  EnumParOut) 

REG  Enumeration  EnumParln; 

REG  Enumeration  ‘EnumParOut; 

•EnumParOut  «  EnumParln; 
if  ( !  Func3< EnumParln)  ) 

•EnumParOut  *  Ident4; 
switch  (EnumParln) 

{ 

case  Identl:  ‘EnumParOut  =  Identl;  break; 

case  Ident2 :  if  (IntGlob  >  100)  ‘EnumParOut  *  Identl; 

else  ‘EnumParOut  -  Ident4; 
break ; 

case  Ident3:  ‘EnumParOut  *  Ident2;  break; 

case  Ident4:  break; 

case  Ident5;  ‘EnumParOut  -  Ident3; 

} 

} 

Proc7( IntParll,  IntParI2,  IntParOut) 

OneToFifty  IntParll; 

OneToFifty  IntParI2; 

OneToFifty  ‘IntParOut; 

REG  OneToFifty  IntLoc; 

IntLoc  •  IntParll  +  2; 

‘IntParOut  »  IntParI2  +  IntLoc; 


Proc8(ArraylPar,  Array2Par,  IntParll,  IntParI2) 

ArraylDim  ArraylPar; 

Array 2Dim  Array2Par; 

OneToFifty  IntParll; 

OneToFifty  IntParI2; 

REG  OneToFifty  IntLoc; 

REG  OneToFifty  Intlndex; 

IntLoc  ■  IntParll  +  5; 

ArraylPar [IntLoc]  ■  IntParI2; 

ArraylPar [ IntLoc+1]  *  ArraylPar [ IntLoc] ; 

ArraylPar [ IntLoc+30]  »  IntLoc; 

for  (Intlndex  *  IntLoc;  Intlndex  <»  (IntLoc+1);  ++Intlndex) 
Array2Par [IntLoc] [Intlndex]  -  IntLoc; 

++Ar ray2Par [IntLoc] [ IntLoc-1] ; 

Ar ray2Par [ IntLoc+20] [ IntLoc]  »  ArraylPar [ IntLoc] ; 

IntGlob  *  5; 
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Dhrys  tone 


Enumeration  Fund ( CharPar 1 ,  CharPar2) 

Capi talLetter  CharParl; 

Capi talLetter  CharPar2; 

REG  Capi talLetter  CharLocl; 

REG  CapitalLetter  CharLoc2; 


} 


CharLocl  »  CharParl; 
CharLoc2  *  CharLocl; 
if  (CharLoc2  !*  CharPar2) 
return  ( Identl ) ; 

else 

return  ( Ident2  )  ; 


boolean  Func2 ( StrPar I 1 ,  StrParI2) 

String30  StrParll; 

String30  StrParI2; 

REG  OneToThirty  IntLoc; 

REG  CapitalLetter  CharLoc; 

IntLoc  «  1; 
while  (IntLoc  <«  1) 

if  ( Fund ( StrPar 1 1 1 IntLoc] ,  StrPar 12 [ IntLoc  +  1] ) 

CharLoc  ■  'A'; 

♦♦IntLoc ; 

} 

if  (CharLoc  >■  'W  tb  CharLoc  <*  'Z') 

IntLoc  »  7; 
if  (CharLoc  'X') 

return (TRUE) ; 

else 

{ 

if  ( strcmp(StrParIl,  StrParI2)  >  0) 

IntLoc  ♦■  7; 
return  (TRUE); 

} 

else 

return  (FALSE); 


boolean  Func3(EnumPar In) 

REG  Enumeration  EnumParln; 

{ 

REG  Enumeration  EnumLoc; 

EnumLoc  *  EnumParln; 

if  (EnumLoc  »»  Ident3)  return  (TRUE); 
return  ( FALSE) ; 


*  Identl 
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Dhry stone 


lifdef  MOSTRUCTASSIGN 
memcpy ( d ,  s,  1) 


register  char 
register  char 
register  int 

while 

) 

lend  i  f 


*d; 

*«; 

1; 

(1--)  »d++  - 


*SH; 


4-4 


Block  Write 


/* 

*  bloc  kwr 1 te . c 

* 

*  This  program  creates  a  very  large  file. 

* 

*  Compile  by:  cc  -O  blockvrite.c  -o  blockvnte 

V 


•define  NAME 
•define  FNAME 
•define  BSIZE 
•define  BLOCKS 


"BLOCKWRITE" 

"bigfile" 

8096  /*  8K  block  */ 

128  /*  number  of  8K  blocks  in  file  ( tota.  1  Mbyte)  */ 


•include  <sys/file.h> 


i n  t  main()  { 

int  f i lef lags  *  OCREAT  OTRUNC  0_APPEND  0_WRONLY; 

int  filemode  *  0777; 

i  n  t  f  ; 

int  lcoun t  »  0 ; 

char  buffer [BSIZE] ; 

int  i  ; 

prmtf("%s:  beginning  ( %d  bytes  in  file))tn”,  NAME,  BSIZE  *  BLOCKS) 
f flush! 1); 

for(lcount  *  0;  lcount  <  23;  lcount++)  ( 

if  ((f  »  open ( FNAME ,  fileflags,  filemode))  <»  0)  { 

printf("%s:  unable  to  create  '%s'»n*,  NAME,  FNAME) ; 
exit(l); 

} 

for  ( i-1; i< -BLOCKS; i++)  write(f,  buffer,  BSIZE); 
close! f ) ; 

} 

printf(*%s:  complete  ( %d  bytes  in  file))*n",  NAME,  BSIZE  *  BLOCKS); 
f flush! 1) ; 


Block  Read 

/* 

*  olockread  c 

* 

*  This  proqram  reads  a  very  large  file, 

*  Comp)ie  by-  cc  O  ftiorkread.c  o  biockread 

a 

*  / 

♦  de  f ;ne  NAME  ’BLOCKREAD" 

♦def.ne  ENAME  "bigfiie" 

♦define  BSIZE  8096  /*  8K  bloc*  •/ 

♦define  BLOCKS  128  /•  number  of  8K  blocks  in  file  (total  1  Mbyte)  */ 

♦  include  •  sys/f  ;*.e  n> 

int  ma ; n (  )  { 

int  file  flags  =  0_RDON  L i  ; 

int  f liemode  ■  04  44  ; 

int  f  , 

char  but  f  er  (  BS  I  7.E)  ; 

;  n  r 

int  ICOUi'f  ; 

pr i nt  f  (  *  %* .  beginning  (  %d  bytes  in  fiie))»n*,  NAME,  BSIZE  *  BLOCKS) 
f  f lu«h< 1 ) ; 

torjlcount  ■  0;  lcount  <  39;  lcount++)  { 

if  ( ( f  *  open { FRAME ,  fileflags,  filemoda))  <»  0)  { 

printf("%s:  unable  to  open  '%«'*n",  NAME,  FNAME ) ■ 
e*  i  t ( 1 )  ; 


for  ( i « 1 ; i < -BLOCKS ; i-M- )  read(f,  buffer,  BSIZE); 
c lose ( f ) ; 

\ 

printfi"%s;  complete  ( %d  bytes  in  file))«n",  NAME,  BSIZE  *  BLOCKS); 

f  f  i  u  s  h  N  ) ; 


Sort  (Part  1  ) 


tcsh  script  to  run  timing  on  sort 

I 

t  sorttest 

I  sortfile  -  input  to  sort 

I  sortstandard  -  presorted  output  for  checking 

» 

echo  Start  of  sort 

sort  -4  + 5  . /preavard/sor t f i le  >  sortout 
echo  End  of  sort.  Start  of  compare, 
diff  . /preaward/sor ts tandard  sortout 
echo  End  of  compare, 
rm  - f  sor  tout 
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Sort  (Part  11) 


.  de  sh 
.  br 
.  ne  5 
.  PP 

• fB« •$ 1 • f R 
PP 


if 

n 

.  ds 

ua 

•  o  ' 

if 

t 

.  ds 

ua 

c  (  aa 

if 

n 

.  ds 

aa 

• 

if 

t 

.  ds 

aa 

•  ( aa 

if 

n 

.  ds 

ga 

1  f 

t 

.  ds 

ga 

i  ( ga 

if 

► 

.  tr 

*•( 

*  * 

.TH  CSH  1  "18  July  1983" 

.  UC  4 

SB  NAME 

csh  •-  a  shell  (command  interpreter)  with  C-like  syntax 
.SH  SYNOPSIS 
.  B  csh 
[ 

8  i-cef •  instvVxX 

1  ( 

arq  . 

1 

.SH  DESCRIPTION 
.  I  Csh 

is  a  first  implementation  of  a  command  language  interpreter 
incorporating  a  history  mechanism  (see 
B  "History  Substitutions)" 

]ob  control  facilities  (see 
.B  Jobs) 

and  a  C-like  syntax. 

So  as  to  be  able  to  use  its  job  control  facilities,  users  of 
.  I  csh 

must  (and  automatically)  use  the  new  tty  driver  fully  described  in 
, IR  tty  ( 4 ) . 

Th ; s  new  tty  driver  allows  generation  of  interrupt  characters 
from  the  keyboard  to  tell  jobs  to  stop,  See 
.  IR  st  ty  (  I ) 

for  details  on  setting  options  in  the  new  tty  driver. 

.  PP 

An  instance  of 
.  I  csh 

begins  by  executing  commands  from  the  file  .cshrc'  in  the 
. I  home 

directory  of  the  invoker. 

If  this  is  a  login  shell  then  it  also  executes  commands  from  the  file 
. login '  there . 

It  is  typical  for  users  on  crt's  to  put  the  command  stty  crt '  '  in  the 
.1  ifc  . login 

file,  and  to  also  invoke 
IR  tset.  (  l ) 
there . 

.  PP 


Sort  (Part  11) 


In  the  normal  case,  the  shell  will  then  begin  reading  commands  from  the. 
terminal,  prompting  with  % 

Processing  of  arguments  and  the  use  of  the  shell  to  process  files 
containing  command  scripts  will  be  described  later. 

.PP 

The  shell  then  repeatedly  performs  the  following  actions: 
a  line  of  command  input  is  read  and  broken  into 
. IR  words  . 

This  sequence  of  words  is  placed  on  the  command  history  list  and  then  parsed. 
Finally  each  command  in  the  current  line  is  executed. 

.  PP 

When  a  login  shell  terminates  it  executes  commands  from  the  file  .logout' 
in  the  users  home  directory 
. sh  "Lexical  structure" 

The  shell  splits  input  lines  into  words  at  blanks  and  tabs  with  the 
following  exceptions. 

The  characters 
i  ‘  '  ;  '  < ’  > '  ('  )' 

form  separate  words. 

If  doubled  in  it',  «ii  << '  or  >>’  these  pairs  form  single  words. 

These  parser  metacharacters  may  be  made  part  of  other  words,  or  prevented  their 
special  meaning,  by  preceding  them  with  '  se 1  . 

A  newline  preceded  by  a  ie’  is  equivalent  to  a  blank. 

.  PP 

In  addition  strings  enclosed  in  matched  pairs  of  quotations, 

•* ( aa ’ ,  «* (ga ’  or  " ’ , 

form  parts  of  a  word;  metacharacters  in  these  strings,  including  blanks 
and  tabs,  do  not  form  separate  words. 

These  quotations  have  semantics  to  be  described  subsequently. 


Sort  (Part  111) 


$  last  argument 

&  Repeat  the  previous  substitution. 

0  first  (command)  word 

10  es  wr i te . c 

11  cat  oldwrite.c 

12  d  i  f f  *wr i te , c 

[1!  1234 

•*(ua  first  argument,  i.e.  1' 

•-•flyifR  abbreviates  0 » - » f Iy • f R« ! 1 

•09  write  michael 

•flntfR  «fIn*fR»  1 th  argument 

•flxifR*1*  abbreviates  iflxtfRi  •-$' 

• f Ixi f R*  • - • f Iy • f R  range  of  words 

d  directory 

e  existence 

f  plain  file 

o  ownership 

r  read  access 

s/« f I 1 • f Ri i /• f Ir • f R* I /  Substitute  sfllifR  for  iflrifR 
w  write  access 

x  execute  access 

z  zero  size 

$  f  name 

$$ 

$* 

50 

$< 

$?0 

S  ?name 
Sname 

Sname [ selector ] 

Snumber 
$  f | name  > 

${ ?name) 

${ name [ selector ]  } 

$ ( name } 

S t  number  } 

(4s  in 
( So  th 

(See  the  description  of 
(The 

(The  words 
(as  in 

(e . g .  ' $shell '  ) . 

( second  form) . 

-B  «-V 
.B  i-X 
.B  i-c 
.  B  « -e 

.8  »-f 

.  B  •  -  i 
.  B  •  -  n 
.8  •  -s 

.B  *-t 


Sort  (Part  111) 


.B  l-v 
•  B  i-x 
•B  alias 

.B  alloc 
.B  break 
.B  breaksw 
.B  breaksw 
.B  breaksw 
.B  cd 
.  B  chdir 
.B  continue 
.B  default: 
B  default: 
B  else 
B  else 
B  end 
B  end 
B  end 
B  end  if 
B  end  if 
B  endsw 
B  endsw 
B  exit 
B  history 
B  login 
B  logout 
B  nice 


Integer  Arithmetic 


/* 

*  integer  .  c 

* 

*  This  program  does  integer  arithmetic 

* 

1  Compile  by:  cc  -O  integer, c  -o  integer 

* 

*/ 


♦define  NAME  "INTEGER" 

♦define  COUNT  2900000  /*  number  of  iterations  */ 


main(  )  { 

long  i;  /*  iteration  counter  */ 

long  a,  b,  c,  d;  /*  integer  variables  for  arithmetic  */ 


a  *  1234;  b  =  2345;  c  *  3456;  d  «  4567; 

printf(*%s:  beginning  ( %d  i terations ) *n" ,  NAME,  COUNT); 
for  (  i=0;  i<COUNT;  i++)  {  /*  do  some  arithmetic  */ 

a  =  b  c  •  d ; 
b  --  a  *  b  /  d; 

} 

printf("Ss:  complete  ( %d  i terat ions ) «n" ,  NAME,  COUNT); 
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Real  Arithmetic 


/* 

*  real.c 

* 

*  This  program  does  real  arithmetic. 

* 

*  Compile  by:  cc  -O  real.c  -o  real 

* 

*/ 

tdefine  NAME  "REAL" 

♦define  COUNT  600000  /*  number  of  iterations  */ 

float  aa,bb,cc,dd; 

int  i  i  , j  j , k k ; 

main(  ) 

{ 

pr  int  f ( " %  s  :  beginning  ( %d  iterations ) in" ,  NAME,  COUNT) 
forfii  «  1;  ii  <  COUNT;  ii++)  { 
aa  *  i i ; 

bb  «  aa  •  aa; 

cc  *  ( bb  -  aa  -  ,137526)/aa; 
pr intf ( "%s :  complete  (%d  iterations ) in* ,  NAME,  COUNT); 


Large  Data  Space 


/* 

*  largedata.c 

* 

*  This  program  has  a  data  space  Larqer  than  real  memory 

* 

*  Com^ilir  by;  cc  -O  largedata  c  o  largedata 

* 

*/ 


♦define  NAME  "LARGEDATA" 

♦define  COUNT  20000  /*  number  of  iterations  */ 


♦define  BLOCK 
♦define  BSIZE 
♦define  AODR 


1024  /*  IK  block  */ 

4000  /*  large  buffer  sue  (blocks)  •/ 

OxeOOO  /*  base  address  of  array  */ 


main(  )  { 

register  char 
register  long 
register  long 
register  long 
register  long 


•curptr ; 
i; 

pagecount ; 

limit; 

size; 


/*  current  pointer  */ 

/*  iteration  counter  */ 

/*  number  of  new  pages  */ 
/*  number  of  references  */ 
/*  size  of  array  */ 


limit  «  COUNT; 
size  »  BSIZE; 


sbr k ( AOOR  +  BSIZE  *  BLOCK);  /*  increase  data  space  */ 

srand(l);  /*  init  random  generator  */ 

i  »  0; 

pagecount  *  0; 

printf("4s;  beginning  ( %d  iterations,  size  %d)»n", 

NAME,  COUNT,  size  *  BLOCK); 

while  (  *+i  <  limit  )  {  /*  make  COUNT  memory  references  */ 
curptr  *  (char  * ) ( ADDR  +  ((rand()  %  size)  *  BLOCK)); 
if  (  *curptr  *»  0  )  { 

pagecount+i;  /*  increase  new  page  count  */ 

•curptr  *  1; 

> 

printf("%s;  complete  ( %d  pages  referenced,  %d  for  the  first  time) in”, 
NAME,  COUNT,  pagecount); 


Compile 


t  Compile  and  load  of  to  routine 

t 

t  compiletest 
I  to.c  -  C  source 
I  subs.c  -  C  source 

* 

echo  cc  -O  -c  preaward/to . c 
cc  -0  -c  preaward/to . c 
echo  cc  -O  -c  preaward/subs . c 
cc  -0  -c  preaward/subs . c 
echo  cc  -O  -o  to  to.o  subs.o 
cc  -0  -o  to  to.o  subs.o 
rm  - f  to.o  subs . o 


Ss 


